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The production of two jets is the simplest exclusive quantum chromodynamics process in electron- 
positron annihilation. Using this process, we examine the structure of next-to-next-to-leading order 
(NNLO) corrections to jet production observables. We derive a subtraction formalism including double 
real radiation at tree level and single real radiation at one loop. For two-jet production, these sub- 
traction terms coincide with the full matrix elements, thus highlighting the phase space structure of 
the subtraction procedure. We then analytically compute the infrared singularities arising from each 
partonic channel. For the purely virtual (two-parton) NNLO corrections, these take the well known 
form predicted by Catani's infrared factorization formula. We demonstrate that individual terms in 
the infrared factorization formula can be identified with infrared singular terms from three- and four- 
parton final states, leaving only single poles and a contribution from the one-loop soft gluon current, 
which subsequently cancels between the three- and four-parton final states. Summing over all different 
final states, we observe an explicit cancellation of all infrared poles and recover the known two-loop 
correction to the hadronic i?-ratio. 



Abstract 



1 Introduction 



Jet production observables can be measured very accurately at present high energy colliders. By con- 
fronting these data with theoretical calculations, one can determine the strong coupling constant or parton 
distribution functions [1]. Analyzing the different sources of error on these determinations, it becomes 
clear that the largest source of uncertainty are currently the insufficiently precise theoretical calculations 
of QCD corrections to jet observables, which mostly include terms only to the next-to-leading order (NLO) 
in perturbation theory. 

To improve this situation, the calculation of next-to-next-to- leading order (NNLO) corrections to jet 
observables is mandatory. For an n-jet observable, several ingredients are required; the two-loop n-parton 
matrix elements, the one-loop (n+l)-parton matrix elements and the tree level (n+2)-parton matrix el- 
ements. In the recent past, enormous progress has been made in the calculation of two-loop 2 — ► 2 and 

1 — > 3 QCD matrix elements, which are now known for all massless parton-parton scattering processes [2-5] 
relevant to hadron colliders as well as for 7* — > qqg [6,7] and its crossings [8]. For the corresponding par- 
tonic processes, the one-loop matrix elements with one additional parton [9, 10] and the tree-level matrix 
elements with two more partons are also known and form part of NLO programs for 1 — * 4 [11, 12] and 

2 — > 3 reactions [13]. Since these matrix elements lead to infrared singularities due to one or two partons 
becoming theoretically unresolved (soft or collinear), one needs to find one- and two-particle subtraction 
terms which account for these singularities in the matrix elements, and which arc sufficiently simple to 
be integrated analytically over the unresolved phase space. One-particle subtraction at tree level is well 
understood from NLO calculations [14-18] and general algorithms are available for one-particle subtraction 
at one loop [19-22], in a form that has recently been integrated analytically [21,22]. 

Tree-level two-particle subtraction terms have been extensively studied in the literature [23-25] . However 
their integration over the unresolved phase space remains an outstanding issue. The problem of integrating 
out double real emission contributions has so far only been addressed in specific calculations [26-29] , each 
of which requires a subset of the ingredients needed for generic jet observables at NNLO. 

In this paper, we aim to contribute to a formulation of a general subtraction scheme including double 
real emission by studying the simplest QCD process, two-jet production in e + e~ annihilation at NNLO. 
Using the iterated sector decomposition [30,31] for the treatment of all real emission contributions [32-34] 
and without the need to define subtraction terms, the NNLO corrections to this process have been obtained 
in a numerical form [35] very recently. 

In this work, we extract the infrared divergent parts of all partonic contributions to two-jet final states 
analytically and identify them with known structures in the virtual two-loop [36] and one-loop single 
unresolved [20] corrections. These identifications might help to construct a general NNLO subtraction 
formalism, which would allow the NNLO computation of three-jet production in e + e~-annihilation and, 
in the longer term, the inclusion of subtraction terms for double initial state emission, as required for the 
NNLO corrections to jet observables at lepton-hadron and hadron-hadron colliders. 

This paper is organized as follows. In Section^ wc describe the different partonic contributions yielding 
two-jet final states at NNLO. Section describes a formalism to subtract the infrared singularities from 
single and double real emission present in these contributions. All partonic channels are computed in 
Section 01 and decomposed into infrared finite and infrared divergent parts, the latter being identified with 
known structures. The integrated partonic terms represent the analytic integrals of the subtraction terms 
accounting for all real singularities. The structure of the infrared cancellations is illustrated in Section [5] 
where we also check the correctness of our results by rederiving the NNLO corrections to the hadronic 
i?-ratio. An appendix listing all integrals needed in this calculation is enclosed. 



2 Perturbative corrections to two-jet final states 

Two-jet final states in e + e~ annihilation are produced by the primary process 7* — > qq, the decay of a 
virtual photon into a quark-antiquark pair, 

7*(g)— > g(Pi) + g(pa). (2.1) 
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At higher orders in perturbation theory, this process receives corrections from the exchange of virtual 
particles and from real radiation. While the NLO corrections involve up to three final state partons, one 
finds at NNLO that two-jet final states are produced through partonic final states including up to four 
partons. The individual partonic channels are: 
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The contributions of these processes are weighted by so-called jet functions (described in detail in 
Section below), which select two-jet final states from the partonic final state momenta. For the ana- 
lytic extraction of infrared poles present in all partonic channels, we will determine the inclusive z-parton 
production cross sections for these processes. They are described by a single Lorentz-invariant 




The n-particle phase space in dimensional regularization with d = 4 — 2e space-time dimensions reads, 

d$ n = P \ - . . . Q P " (27r) d S d (q - Pl -...-p n ). (2.3) 

Parameterizations of the three- and four-particle massless phase space appropriate to analytic integration 
are discussed in [33] . In the following, we frequently normalize the three- and four-particle phase space to 
the volume of the two-particle phase space, 

P 2 = j d$ 2 = 2-^ n-^ feV ■ (2-4) 



3 Infrared subtraction terms 

To obtain the perturbative corrections to a jet observable at a given order, all partonic multiplicity channels 
contributing to this order have to be summed. Each partonic channel contains infrared singularities. 
However, after summation, all singularities cancel among each other [37]. 

While infrared singularities from purely virtual corrections are obtained immediately after integration 
over the loop momenta, their extraction is more involved for real emission (or mixed real-virtual) contri- 
butions. Here, the infrared singularities become only explicit after integrating the real radiation matrix 
elements over the phase space appropriate to the jet observable under consideration. In general, this in- 
tegration involves the (often iterative) definition of the jet observable, such that an analytic integration is 
not feasible (and also not appropriate) . Instead, one would like to have a flexible method that can be easily 
adapted to different jet observables. Therefore, the infrared singularities of the real radiation contributions 
should be extracted using infrared subtraction terms. These subtraction terms are constructed such that 
they approximate the full real radiation matrix elements in all singular limits while still being sufficiently 
simple to be integrated analytically over a section of phase space that encompasses all regions corresponding 
to singular configurations. 

To specify the notation, we define the tree level m-parton contribution to the J-jet cross section in 
c?-dimcnsions by, 

da B = N m J2^ m (pi,---,p m ,Q)^\M m (pi,... Pm )\ 2 ^j m \pi,... lPm ). (3.1) 
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A/i n includes all QCD-independent factors, denotes the sum over all configurations with to partons, 
d<fr m is the phase space for to partons, S m is a symmetry factor for identical partons in the final state and 
finally |A4 m | is the tree level m-parton matrix element. The jet function T^J"^ defines the procedure for 
building J-jets out of to partons. The main property of Jj is that the jet observable defined above is 
collinear and infrared safe as explained in [18]. In general contains 9 and (5-functions. and J 7 ^ = 1. 



3.1 NLO infrared subtraction terms 

At NLO, we consider the following m-jet cross section, 



da nlo 



d* m+ i 



(dc 



NLO 



d °NLo) 



dcr NLO + 



d*„ 



d$„ 



da NLO 



(3.2) 



The cross section dc m 
m — ► m + 1 



while da NLO 



NLO has the same expression as the Born cross section da NLO (|3.1|l above except that 



cross section dcr^ LO 



as da^Lo m a ^ appropriate limits 



is a local counter-term for do~fj L0 . 



is the one-loop virtual correction to the m-parton Born cross section da B . The 

It has the same unintegrated singular behavior 
Their difference is free of divergences and can be integrated over the 
(to + l)-parton phase space numerically. The subtraction term d<7^ iC) has to be integrated analytically 
over all singular regions of the to + 1-parton phase space. The resulting cross section added to the virtual 
contribution yields an infrared finite result. 

The subtraction term da^LO can be constructed in a number of different ways, and even the phase space 
used with these subtraction terms can vary from method to method. In general, two different methods 
are applied at next-to-leading order. In the slicing method [15, 17], one exploits the fact that both real 
radiation matrix element and final state phase space factorize in all soft and collinear limits. Therefore, a 
subtraction term for a given singular limit can be obtained by simply expanding the full matrix clement 
around the limit under consideration. This subtraction term is then integrated over a small slice of phase 
space, which is tailored to contain only one singular configuration. In the subtraction method [14, 16, 18], 
the infrared subtraction terms are integrated over the whole final state phase space (or at least over the 
full sub-phase space appropriate to one or more final state momenta). In this method, care has to be 
taken in order to construct subtraction terms which fully account for the limit they are aimed at without 
introducing spurious infrared singularities in other limits. 

A systematic procedure for finding NLO infrared subtraction terms in the second method is the dipole 
formalism derived by Catani and Seymour [18]. Their subtraction terms are obtained as sum of dipolcs 
^2T>ijk (where each dipole corresponds to a single infrared singular configuration) such that, 



da NLO " d(7 NLO 



= N in d$ m+ i(pi,. 



L,Q)- 



1 



\M m+ l(Pl,-Pm+l)\ 2 ^ r m+1 \ Pl , 



E 

pairs ij k^i,j 



^2 V ijk \M m ((?pl, -Pij,Pk, ■■■,P m +l)\ 2 Tjipi, ~Pij,Pk, -,Pm+l) 



(3.3) 



The dipole contribution T>ijk involves the TO-parton amplitude depending only on the redefined on-shell 
momenta pi, ..pij,pk,p m +i and the splitting matrix Vy,fc/ 'sy which depends only onpi,pj,pk- The momenta 
Pi, pj and Pj. are respectively the emitter, unresolved parton and the spectator momenta corresponding to 
a single dipole term. The redefined on-shell momenta pij,pk are linear combinations of them. 

For the NLO two-jet cross section, the real corrections involve three partons in the final state. The 
single unresolved configurations occur when the invariants S13 or S23 vanish. We notice that the subtraction 
term built as the sum of the dipoles is very close to the full three-parton matrix element. This argument 
can in fact be turned around: the dipole subtraction terms can be obtained by an appropriate partial 
fractioning of the full matrix element squared into terms which can be uniquely attributed to individual 
single unresolved configurations. In the following, we shall use this formulation as definition of a dipole 
term, without specifying its precise form. 
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The jet function !Fj"^ in l|3.3[l depends not on the individual momenta pi , pj and pu , but only on jjy , pk ■ 
One can therefore carry out the integration over the unresolved dipole phase space appropriate to pi, pj 
and pk analytically, exploiting the dipole factorization of the phase space [18], 

d$ m+ i(pi, ...,p m+ i,Q) = d$ m (pi, ..pij,pk, —,Pm+i) ■ d$ D (pi,pj,p k ) , (3.4) 

where d$D is symmetric under permutations of i, j, k. The dipole phase space is proportional to the three- 
particle phase space, as can be seen by using m = 2 in the above formula and exploiting the fact that the 
two-particle phase space is a constant, 

p 2 = y"d$ 2 , (3.5) 

such that 

d$ 3 = P 2 d$ D . (3.6) 

It should be noted in this context, that the sum of dipoles is much easier to integrate over the dipole 
phase space than an individual dipole. In what follows, we shall take into account this fact, and choose 
to use the full three-parton matrix element as subtraction term [12] to the real corrections to the two-jet 
rate at NLO (and, in Section l3~2l below, the full four-parton matrix element as candidate subtraction term 
at NNLO). Taking the full matrix element as the subtraction term provides a clear visualization of the 
phase space regions covered by the different contributions involving the jet-functions . At NLO the 
difference of the unintegrated real and subtracted contributions to the two-jet cross section reads, 

dcr NLO ~ dcr NLO = N ln d<P 3 { Pl ,...,p 3 ,Q)\M 3 {p l7 ...p 3 )\ 2 x 



(3.7) 



(21 

The factors 1/2 in front of the T\ are obtained by first writing out both dipoles, i.e. both terms obtained 
by partial fractioning of the three-parton matrix element which are appropriate to the configurations where 

(2) 

parton 3 is unresolved with respect to either parton 1 or 2. Exploiting the fact that !F§ ' = 1, these dipoles 
can be added together to yield the full matrix element. Reintroducing both jet functions consequently 
yields a factor 1/2 for each of the momentum redefinitions. 

Equation (|3.7|l can be interpreted as follows: out of the three-parton final state phase space one considers 

(3) (2) 

the two-jet configuration (defined by T 2 ) minus the fully inclusive phase space (defined by T\ — 1), 
which covers to this order the two- and three-jet configurations. 
For the analytic integration, the subtraction term is rewritten as, 

da s NLO = d$ 3 \M 3 \ 2 = d<S> 2 \M 2 \ 2 I \M 3 \ 2 , (3.8) 

J d<i>D 

where the three-parton phase space is now a product of the dipole phase space and a two-parton phase 
space and where the three-parton matrix element is normalized to the two-parton matrix element such that 

The analytic integral of the subtraction term is therefore the three-parton contribution integrated over the 
fully inclusive phase space. The three-parton contribution to two-jet final states is thus obtained by adding 
and subtracting the three-parton contribution to the inclusive cross section. The subtracted term is then 
used in the numerical integration, accounting for all infrared singularities in the two-jet final states, while 
the added term is integrated analytically to make the infrared singularities explicit. 

To summarize schematically the contributions to the two-jet cross section at NLO, one finds the following 
two- and three-parton phase space contributions 

da 2 N i =d$ 3 \M 3 ( Pl ,... P3 )\ 2 \(4 3) -4 2) )} +^2\M 2 \ 2 T (2) ([ |M 3 | 2 + |MP| 2 Y (3.10) 

where | Af^' 1 1 2 is the normalized one- loop virtual correction to the j-parton matrix element, which exactly 
cancels all infrared poles obtained from the integral of the subtraction term. Both virtual and subtraction 
terms are proportional to the two-parton matrix element and phase space. 
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3.2 NNLO infrared subtraction terms 



At NNLO, the two-jet production is induced by final states containing up to four partons, including the 
one-loop virtual corrections to three-parton final states. As at NLO, one has to introduce subtraction terms 
for the three- and four-parton contributions. Schematically the NNLO two-jet cross section reads, 



d^NNLO - / ( d(J NNLO ~ d(7 NNLo) + / da NNL( > 



d(T NNLO ~ d<T NNLo) + / d<T NNLO + d<T NNLO > ( 3 ' U ) 

/d<I>3 N ' Jd$3 Jd<I>2 

where da% NLO denotes the real radiation subtraction term coinciding with the four-parton tree level cross 
section da^ NLO in all singular limits. Likewise, da^^Q is the one-loop virtual subtraction term coinciding 
with the one-loop three-parton cross section da^^ L Q in all singular limits. Finally, the two-loop correction 
to the two-parton cross section is denoted by da^ NL0 . 

A method to construct subtraction terms for n-parton cross sections at NNLO has been worked out in 
detail by Weinzierl for double real radiation [25] and single real radiation off one- loop matrix elements [22]. 
The method which we will outline in the following differs from that proposed in [22,25], and we shall 
highlight the main differences at the end of this section. 

Following the line of thought elaborated in the NLO discussion above, the subtraction terms are chosen 

to be the full matrix elements, i.e. the sum of the full four-parton matrix elements related to the tree 

level processes 7* — > qqgg, 7* — > qqqq and 7* — > qqq'q'(q 7^ q') for da s and the one-loop corrected matrix 

v s 1 

element squared of 7* — > qqg for d<r N ^ L0 . 

Clearly, these subtraction terms have the appropriate behavior in all singular limits. Moreover, as 
explained in [33] , they can be integrated analytically by reducing them algebraically to four master integrals 
which have been derived explicitly in [33] . Likewise the integration of the one-loop correction to 7* — > qqg 
over the dipolc phase space can be performed analytically by reducing it to three master integrals which 
are listed in the appendix. These calculations are presented in Section 0] below. 

To perform the phase space integration of the subtraction terms, one needs to factorize the phase 
space into a product of a two-particle phase space multiplying a dipole phase space for the three-particle 
contribution and multiplying a tripole phase space 

d$ 4 = Pi d^T, (3.12) 

for the four-particle contribution. The latter factorization [24, 27, 33] is obtained by redefining a set of 
four massless on-shell momenta (emitter, two unresolved partons, spectator) into two on-shell massless 
momenta. 

The construction of the NNLO subtraction term to the four-parton tree level cross section starts from 
the full four-parton matrix element: 

da S / NLO =N in d^(p 1 ,... > p 4 ,Q) \M A ( Pl ,... P4 )\ 2 Q^ 2 (2) (pi34,p 2 ) + i^ 2) (p 23 4,Pi)) , (3.13) 

which is in analogy to the NLO subtraction term. However, in order to have a subtraction term which 
correctly accounts for all singularities in the two-jet region and can be integrated over the inclusive (i.e. 
two-jet, three-jet and four-jet) phase space, this is not sufficient. The infrared singularities yielded by the 
subtraction term in the two-jet region of the four-parton phase space exactly cancel those coming from the 
double real emission matrix element integrated over the same two-jet region. Furthermore, the subtraction 
term yields only finite terms when integrated over the four-jet region. However, the subtraction term 
generates single infrared singularities when integrated over the three-jet region of the four-parton phase 
space which need to be treated carefully. Indeed, the singular structure in the three-jet region has to be 
mapped out and subtracted from the subtraction term. This is done by subtracting the sum of all dipole 
terms appropriate to single-particle singularities from the four-parton matrix element subtraction term. 
According to (|3.3(l . these read 

da NNLO = N ind®4,(pi,—,P4,Q) J! ^ V ijk\M 3 (j)l,Pij,Pk)\ 2 ^3 (pi,Pij,Pk) , (3-14) 

pairs i,j k^i,j 
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where I k. With this, the full NNLO subtraction term becomes 

S SO 5 1 

dCjViVLO = ^ a NNLO ~ d<J NNLO • 

Together with 

d^NNLO = N in d$ 4 (Pl , ■ ■• , Pi, Q) \ MiiPl , .. -Pi) | 2 ^ 4) (pi , p 2 , P3, Pi) 

we obtain the difference between double real and subtracted matrix element as, 

N in d$ 4 (pi, -,Pi,Q) x 



(3.15) 
(3.16) 



aa NNLO 



da S 

aa NNLO 



M \(p\. ...p;) 1 [ ^ 4) (pi,p 2 ,P3,P4) - ^•? 7 2 2) (Pl34,P2) - TjflP (P234, Pi) 



,Pij,Pk) 

pairs i,j k^i,j 



(3.17) 



which denotes the full contribution from the four-parton channel to the two-jet production cross section. 
It is finite and can therefore be integrated numerically. 

The subtraction term for the one-loop virtual contribution is constructed in complete analogy to the 
NLO three-particle subtraction term l|3.7|l by using the full one-loop three-particle matrix element. We 
have, 



da 



NNLO 



da 



VS.l 

NNLO 



Ni„d9!3(pi,...,P3,Q) IM^ipu.-.pa)^ x 



7^2 2) (Pl3,P2) ~ ^2 2) (P23,Pl) 



(3.18) 



The difference between the one- loop virtual contribution and its subtraction term, da^ NLO — &&nnloi is 
integrated over the two-jet and three-jet regions of the three-parton final state phase space. In the two-jet 
region, it is finite (in fact, yields a zero contribution) and can be integrated numerically in a straightforward 
manner. Although dcrj^^ will not produce infrared poles from real emission in the three-jet region of 
the final state phase space, it will still generate infrared poles from this region, since \M.%' 1 {pi, ...ps)\ 2 
contains explicit infrared poles from the loop integration. These poles are proportional to the tree level 
three-particle matrix element squared, and are canceled by extracting the infrared poles from the real 



emission subtraction term da 



S,l 

NNLO' 



To see this, we express 



da 



s.i 

NNLO 



= AT in d$ 4 (pi,...,p4,Q) Y E Vi * k l^iPhPijjPkjl 2 ^a(j>i,Pij,Pk) 

pairs i,j k^i,j 



N.„ 



E E d ^3(j>i,Pij,Pk,Q) \M a (pi,pij,p k )\ 2 A f'(pi,P~ij,Pk) / f». 

,irs i 4 h=Li ■> Jd^D 



k .(3.19) 



pairs i,j k^i,j 



Since the sum of integrals of the dipole terms over the dipole phase space exactly cancels the infrared 
singularities present in the corresponding one-loop matrix element, we find that 



da 



VS,1 
NNLO 



da 



S,l 

NNLO 



N in d$ 3 {p 1 ,p 2 ,P3,Q) X 



\M^(p U -P 3 W 



-^2\Pl3,p2) ~ ^2 2 \P23,Pl) 



\M 3 (pi 1 P2 1 P3)\ Z FP(P1,P2,P3) E 

pairs i,j;k^i,j 



v 



d<S> D 



(3.20) 



is finite if integrated over the three-jet region of the three-parton final state phase space. The full contri- 
bution from the three-parton channel to two-jet final states is therefore 



da 



V,l 

NNLO 



da 



VS.l 
NNLO 



da 



S,l 

NNLO 



G 



.1-1 / T@)U. „.\ 1 t(2)^„^^ 1 t( 2 )i 



\m \ \ Pl , . . .p 3 ) r ( *r (pi ,p 2 , p 3 ) - 2 *r (fca , £2 ) - ~n > , w ) 



AT in d$ 3 (pi,p 2 ,P3,Q) 

-|-M 3 (Pl,P2,P3)| 2 ^3 3) (Pl>P2,P3) X! / V% i k 

pairs i,j;k^i,j ^ 

which is finite over the full three-parton phase space and can be integrated numerically. 



(3.21) 



Finally, the subtraction terms da 



S,o 

NNLO 



and da 



vs,i 

NNLO 



integrated over the tripole and dipole phase space 



cancel the infrared singularities present in the two- loop virtual contribution da^ NL0 . Using T { 2 ] = 1, we 
find 



^nnlo = 4d$ 4 \M 4 \ 2 = N m d<P 2 \M 2 \ 2 I |M 

Jd* T 

^nnIo = N in d$4 \M^\ 2 = iV m d$ 2 |X 2 | 2 f \M. 
Adding these to the two-loop virtual contributions, 



V,li2 
3 



j y,2 , s,o , , vs,i 

^ NNLO aa NNLO + aa NNLO 



N in d<S> 2 \M 2 \ 



\M, 



V,2i2 



d$ 7 



|M 4 | 



|M 



V,l|2 



d$r 



(3.22) 
(3.23) 

,(3.24) 



we find an infrared finite result for the two-parton channel. 

To summarize schematically, the NNLO corrections to the two-jet rate take the following structure 



da 



NNLO 



acr NNLO 



da 



S,0 

NNLO 



da 



S,l 

NNLO 



S,l 



A ' A ' A 

aa NNLO — aa NNLO ~ ® NNLO 

, V,2 , S,0 , t VS,1 1 

aa NNLO aa NNLO ' aa NNLO\ 



(3.25) 



d$ 4 



-d$ 3 



|M 4 



r (4) _ ^.(2) 
•'2 •'2 



£l^ 3 | 2 



D 



d$i 



(3) 



-d$ 2 |Ma| 



"K ,a l a + 


/ |Af 4 | 2 + 


/ |M 3 ^| 2 " 




Jd$ T 


Jd<s> D 



r. 



(2) 



(3.26) 



The two-, three- and four-parton final state contributions are made explicit in this formula as those contain 
respectively the corresponding parton phase space d$i. These correspond to (|3.24l) . H3.21fl and (|3.17|l . 
which are all separately finite and can be integrated numerically over the two-, three- and four-parton final 
space respectively. 

The interpretation of the individual terms is as follows: the infrared singularities of the four-parton 

contribution da^ NLO to two-jet final states are subtracted through the full matrix element, which is 

s o 

integrated over the inclusive phase space, da^ NLQ . Since this term also contains infrared singularities in 

S 1 

the three-jet region, a further subtraction term da^ NLO applicable only in this region, has to be subtracted 

S,o 



off da 



NNLO- 



Infrared singularities from the one-loop virtual correction to the three-parton contribution 



to two-jet final states da^^ LO are subtracted through the full one-loop matrix element integrated over 

VS 1 

the inclusive phase space, da N} ^ L Q. This term yields explicit infrared poles in the three-jet region, which 

s 1 

cancel with the infrared singularities obtained by integrating da^ NLO over the single unresolved phase space 
appropriate to three-jet final states. Finally, all infrared poles produced through the two-jet contribution 
of the inclusive phase space integration of the subtraction terms da 



poles from the two-loop virtual corrections, da 



V.2 

NNLO- 



S,0 
NNLO 



and da 



VS,1 

NNLO 



cancel with the 
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The construction of the subtraction terms presented here differs in one essential point from the proce- 
dure proposed in [25] : we first subtract all double and single unresolved singularities from the real radiation 
contribution da^ NLO using dafj^ LO and afterwards add in the oversubtracted single unresolved singular- 
ities through da^NLQ. In contrast, [25] first subtracts all single unresolved singularities from the real 
radiation contribution da^ NLO and then constructs subtraction terms which account only for the double 
real emission contributions in the (already singly subtracted) matrix element. The double real emission 
subtraction terms constructed this way are more involved than the corresponding terms employed here. In 
particular, the subtraction terms obtained in [25] can not be related to four-particle matrix elements in a 
straightforward manner, which renders their analytic phase space integration more difficult. As a matter 
of fact, the subtraction terms constructed here can be integrated using the methods derived in [33], while 
the subtraction terms of [25] first require the computation of as yet unknown classes of integrals [38] . 

On the other hand, let us note that taking the subtraction term to be the full four-parton matrix element 
squared we clearly do not separate singularities arising from individual partonic configurations. Therefore, 
more work is needed for extending the method derived here to final states with more than two jets, while 
the subtraction terms of [25] readily generalize to final states with higher jet multiplicity. 

In the next section, we analytically compute the infrared poles arising in each parton level contribution 
discussed here. For simplicity, we set Fj' = 1 in the following (corresponding to computing the inclusive 
cross section) and discard the single-particle subtraction term daj^l^ LO , whose infrared structure is well 
known from NLO calculations [18] and whose contributions exactly cancel. 



4 NNLO contributions 

As outlined in Section above, two-jet final states at NNLO accuracy are obtained from the production of 
up to four final state partons. In this section, we present the contributions from two, three and four-parton 
final states, elaborating the structure of infrared singularities for each channel. The integrated partonic 
terms derived in this section represent the integrals of the subtraction terms presented in Section 3 and 
account for all singularities arising from real radiation. 



4.1 Two-parton final states 

The calculation of the virtual two-loop corrections to 7* — + qq (also called quark form factor) in dimensional 
regularization with d = 4 — 2e space-time dimensions was performed long ago [39]. In the following, we 
only list the results, and illustrate the application of the infrared factorization formula. Renormalization 
of ultraviolet divergences is performed in the MS scheme by replacing the bare coupling ao with the 
renormalized coupling a s = a s (/i 2 ), evaluated at the renormalization scale /i 2 



2c 



1 



llN-2N f fa 



Gc 



(4.1) 



where 

S f = (47r) e e" e7 with Euler constant 7 = 0.5772 . . . 

and (Aq is the mass parameter introduced in dimensional regularization to maintain a dimensionless coupling 
in the bare QCD Lagrangian density. 

The renormalized amplitude can be written as 



\M) q9 = V4 



7rae 



\M 



(oh 



+ (£) i^ (1) >«*-+(£)V 2 v-+0(ca 



(4.2) 



where a denotes the electromagnetic coupling constant, e q the quark charge, a s the QCD coupling constant 
at the renormalization scale (i, and the are the i-loop contributions to the renormalized amplitude. 

They are scalars in colour space. 

The squared amplitude, summed over spins, colours and quark flavours, is denoted by 



{M\M\ 



\Mh*^qq)\ 2 =T q ,tf). 



(4.3) 
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The perturbative expansion of T qq (q 2 ) at renormalization scale ji 2 = q 2 reads: 



47ra^e 2 



+ 



2tt 



2tt 



T q f(q 2 )+0(al(q 2 )) 



where 



r q f{q 2 ) 



Iqq 
)qq 

-(2)/ 2> 



e)q\ 



{M {0) \M^) qq = AN{\ 
(M ( - 0) \M w ) qq + {MW\MM 



(4.4) 



(4.5) 
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7tt 2 7 A , 
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16+ — + 2 C3 



73^ 4 
1440 



e 2 + 0(e 3 ) 



(M^\M^) qq + (M^\M^) qq + (M^\M^) qq , 



(4.6) 
(4.7) 



where TV is the number of colours. 



In the following, we decompose T q q\q 2 ) into the contributions from the interference of two- loop and 



tree diagrams, 



T 9 ( |< [2X01 V) 



(M^\M^) qq + (M^\M^) 
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1 / 65 



4e 3 9e 2 e V 216 24 
as well as the one-loop self-interference, 
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41tt 
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(4.10) 



The pole structure of T q ^ and of both contributions to T q ^ can be expressed using an infrared factor- 



ization formula, obtained from rcsummation [36]. To this end, we write 

T q f(q 2 ) = Vole Sq f°\q 2 )+^mte q f a \q 2 ), 



Voles 



tf)+finite q f j \q 2 ). 



(4.11) 
(4.12) 
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VoleSqq contains infrared singularities that will be analytically canceled by those occurring in radiative 
processes of the same order. Finite is the renormalized remainder, which is finite as e — > 0. For simplicity 
we set the renormalization scale [x 2 = q 2 . 

The infrared factorization formula shows how to predict infrared pole structure of the one-loop and 
two-loop contributions renormalized in the MS scheme in terms of the tree and renormalized one-loop 
amplitudes, \M.^) and jA^ 1 )) respectively. At one-loop it yields 



Voles 



(lxo) _ 



2ft(M i ° ) \I W {e)\M ( - 0) ), 



(4.13) 



while at two loops 



Toles q f a) = 25? 



+ (M (0) \I {1) {e)\M {1) } 

_ e7 £(l-2e) (0o MR \ (M (0)| J (D (2e) | M (0) ) 



r(l-e) K e 
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(4.14) 



and 



Voles^q* 1 ^ 



2(M W \I {1) (e)\M i()) ) - (M W \I {1) He)I^(e)\M W ) 



where the constant K is 
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(4.15) 



(4.16) 



It should be noted that, in this prescription, part of the finite terms in T^f 1 and T^'^ x ^ are accounted 



ii n 
for by the O(e ) expansion of Voles^*^. 

For the process under consideration, there is only one colour structure present at tree level which is 
just the contraction of the quark and antiquark colours i and j, it is simply 8ij. Adding higher loops does 
not introduce additional colour structures, and the amplitudes are therefore scalars. Similarly, the infrared 
singularity operator J^(e) is a 1 x 1 matrix in the colour space. By explicit evaluation of Voles^ 3 , one 
can show that j' 1 ' (e) is only determined up to a finite constant (which we shall label c\ below) that can 
be chosen arbitrarily. In fact, in [18], where V^'{e) was derived for the first time, a constant term arising 
from the finite contributions to the dipole matrix elements was found. This term was omitted in [36]. 
Including c\, I^ie) is given by, 



2r(i 



N 2 



2N 



-J + - + Cl 



U2 



(4.17) 



(since we have set /i 2 = q 2 ): 



512 



(4.18) 



Note that on expanding S12, imaginary parts are generated, the sign of which is fixed by the small imaginary 
part +i0 of q 2 . Other combinations such as (Ai^ll^ 1 ^ (e) are obtained by using the hermitian conjugate 
operator /^'^(e), where the only practical change is that the sign of the imaginary part of S12 is reversed. 
Finally, the last term of Eq. (|4.14Jl that involves (e) produces only a single pole in e and is given 

by, 



(M (0) \H^(e)\M {0) ) 



4eT(l 



-if( 2 )(M (0) |M (0) ) , 



(4.19) 
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where the constant is renormalization-scheme-dependent, but independent of the choice of c\. As with 
the single pole parts of I^\e), the process-dependent can be constructed by counting the number of 
radiating partons present in the event. In our case, there is only a quark-antiquark pair present in the final 
state, so that 

(4.20) 



where in the MS scheme 
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As a result, one finds the finite parts 



Finite 



(q 2 ) 



4Nq 2 



Tinite 



(2x0), 2 



(q 2 ) 



ANq 2 



A 



£1 
2 

81659 5tt 2 
~ "36" 



(4.21) 
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4.2 Three-parton final states 



(8-cr) 2 
64 



The three-particle final state yielding two jets is 7* 
formed from three partons. 

The renormalized amplitude for this process can be written as 



(4.23) 
(4.24) 

qqg in the configuration where only two jets are 



\M) qqg = V^e q V^cTs [\M^) qqg + (g) \M^) qqg + 0(a 2 s ) 



(4.25) 



where \A4^} are the i-loop contributions to the renormalized amplitude. The squared amplitude, summed 
over spins, colours and quark flavours, is denoted by 



(M\M) qqg = ^\M( 1 * ^^a)\ 2 - 

The perturbative expansion of the squared amplitude at renormalization scale /x 2 = q 2 reads: 

'a s (<7 2 ) 



(4.26) 
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(M^\M^) qqg 



((M^\M^) qqg + (M^\M^) qqg )+0(aUq 2 )) 



(4.27) 



(2) 

The squared amplitudes are integrated over the three-parton phase space and multiplied by ' , the 
jet function ensuring that out of three partons a two-jet final state is formed. The two-jet final state phase 
space for these three-parton configurations contains the regions where the gluon is cither collincar or soft, 
and subtraction terms have to be introduced to extract the resulting infrared singularities from the partonic 
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cross sections. As outlined in Section|21 these subtraction terms are taken to be the full squared amplitudes 
associated with 7* — > qqg up to the one-loop order. Those are integrated over the three-parton phase 
space factorized in a dipole phase space factor d<&£> and the integrated two-parton phase space, P%. The 
momenta associated with this dipole phase space are the two dipole momenta (emitter and spectator) and 
one unresolved (collinear or soft gluon) . The infrared poles present at 0(a s ) are made explicit as follows, 



8vr 2 J d<P D (M (a) \M 
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The infrared poles in T^ g cancel all infrared poles present in the one-loop two-particle final state contri- 

i-qq 



bution in IpTT^ : 



T (4) = Vole S {0x0) 
1 qqg 1 ulca qgg 
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At 0(a 2 ), the renormalized matrix element is (M^^M^qqg + (M^\M^) qq - g . To cancel the poles 
present in the single unresolved contributions related to this virtual matrix element, a subtraction term is 
required, which is chosen to be the virtual matrix element itself. The analytic integration of this subtraction 
term over the dipole phase space is obtained by first reducing all terms in the integral to three master 
integrals, solving iteratively integration-by-parts identities as outlined in [40,41] for two-loop integrals and 
in [33] for four-particle phase space integrals. For this step, extensive use of the algebraic programming 
language FORM [42] is made. As a result, one obtains a linear combination of three master integrals, denoted 
by V5 >a , V5 t b and V&, which are listed in the appendix. Inserting these, one obtains 
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The infrared poles in T Q 
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Z, 
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qqg 



cancel part of the infrared poles from the two-particle final state contribution 
,,,, . This cancellation becomes evident if we decompose 
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and identify individual terms of Voles q q g with terms from the infrared factorization formulae l|4.14|l and 
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The first two terms in this formula cancel with terms in l|4.14[) . while the third term cancels with a term in 
(|4.15(l . The last two terms contain the remaining poles, and will be discussed in detail now. 

For the process under consideration, the infrared singularity operator 1^ (|4.17l) contains only a single 
colour structure (TV 2 — 1)/(2N) = Cf, as does the one-loop squared matrix element . By inserting 1^ 

(|4.17|) in 14.36|) , one finds that all 1/e 4 and 1 /e 3 poles of Xj^l which are proportional to the colour structures 
[(iV 2 - 1)/(27V)] 2 = Cf and {N 2 - l)/(2N)N f = C F N f are correctly accounted for. This cancellation is 
independent of the choice of the constant term c\ in The 1/e 2 terms proportional to these colour 

structures in (|4.36l) depend on c\ . They match the corresponding terms in T}gl if we choose 



43 t: 2 

C1 = T~T 



(4.37) 



With this choice, the first three terms of 1)4. 36|) correctly account for all 1/e 4 to 1/e 2 poles proportional 



to [{N 2 - l)/(2iV)] 2 = Cf and {N 2 - \)/{2N)N f = C F N f in T$], and all remaining single pole terms in 



these two colour structures can be attributed to the H v (e) term in 1)4. 36[) . 

Pole terms corresponding to the colour structure N(N 2 — 1)/(2N) = CfCa can be generated only 
through the combination of and /?o/e in (|4.36() . thus yielding at most 1/e 3 poles. In T^g, one does 
however observe 1/e 4 terms in this colour structure, as first pointed out in [43]. These terms have no 
analogue in the two-particle final state contributions, and must thus cancel among three and four-particle 
final states. We observe that these terms are proportional to the one-loop soft gluon current squared matrix 
element derived in [20] (which is a scalar in colour space for this particular process), integrated over the 
dipole phase space. Explicitly, we find, 
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The integrated soft gluon current accounts for all 1/e to 1/e 2 poles proportional to N(N 2 — 1)/(2N) 

r(2) 
- V 



C f Ca in ^qsgi the remaining single poles are again attributed to Hy\e), which finally reads, 
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The expression for Hy q in the MS scheme is given by, 
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To summarize the discussion of the infrared pole terms arising from three-particle final states at this 
order, we have demonstrated that all double and higher poles can either be identified with terms present 
in the infrared factorization formulae for the virtual corrections (provided an appropriate finite constant 
is inserted in the infrared singularity operator 1^) or with the dipole phase space integral of the soft 
gluon current. The single poles are attributed to a constant Hy \ whose origin (like the origin of the 
appearing in the double virtual corrections) is not fully understood at present. After subtraction of the 
poles, one is left with, 
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4.3 Four-parton final states 

Three different tree level processes can yield four-parton final states relevant to two-jet final states: 7* 
qqq'q' {q 7^ q'), 7* — > qqqq and 7* — > qqgg. The amplitude for each process can be written as 

f(0) 



\M) qqij = ^Anae q Ana s [\M^>) qqij +0(a s )\ , (4.45) 

where ij = q'q',qq,gg. 

The squared amplitude, summed over spins, colours and quark flavours, is denoted by 

(M\M) qq:iJ = ]T |M( 7 * - qqij)\ 2 - (4-46) 

The perturbative expansion of the squared amplitude at renormalization scale — q 2 reads: 

2 



(M\M) qqiJ = 4TraJ24^ 4 



^1) {M^\M^) qqi3 +0{al{q 2 )) 



(4.47) 



These amplitudes need to be integrated over the four-parton phase space and in order to give rise 

(2) 

to contributions to the two-jet rate, they need to be multiplied by the corresponding jet-function T\ . 
Indeed, two of the four final state partons are experimentally unresolved. The two-jet final state phase 
space, however also contain regions where one or two particles are theoretically unresolved; the associated 
partonic cross sections contain single and double singularities. Those are canceled analytically through the 
introduction of subtraction terms. The appropriate subtraction terms for the infrared singular structure 
of these real-real corrections are the corresponding matrix elements themselves as explained in Section |21 
above. These are to be integrated over the tripole phase space factor d<i>T <|3.12(l . which is a factorized 
form of the four-parton phase space d<&4. These integrals have been discussed in detail in [33], where it 
was demonstrated that any integral appearing in this context can be expressed as a linear combination of 
four master integrals R4, Rq, Rs a and Rst, which were derived in [33] and are listed in the appendix for 
completeness. As a result, one finds the integrated subtraction terms, 
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The infrared poles from these four-parton final states exactly cancel the infrared poles present in the 
two- and three-particle final states derived above. To see this cancellation, we consider 
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As for the three-parton final states above, we have to fix c\ according to l|4.37|l to obtain cancellation 
of the 1/e 2 poles in the [(JV 2 - 1)/(2JV)] 2 = Cf and (A^ 2 - t)/(2N)N f = C F N f colour structures. It can 
be seen that the qqgg final state also yields a term canceling the contribution from the one-loop soft gluon 
current in the three-parton channel. Any remaining 1/e poles are collected in 
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The expression for H R in the MS scheme reads, 
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The single pole constants from single virtual and double real radiation add up to the single pole constant 
from the double virtual contribution, 



Subtraction of the pole terms leaves 
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5 Structure of infrared cancellations 



In any infrared safe physical observable, such as the two-jet cross section or the inclusive hadronic production 
rate, one observes the cancellation of all infrared singularities present in individual partonic channels after 
adding all channels contributing to the same (infrared safe) final state [37]. At NLO, this cancellation is 
obvious in 

Voles™ + Voles™ =0. (5.1) 

At NNLO, the cancellation reads 



Voles*** + Voles™ + Voles™ + Voles™ = . 



(5.2) 



In the computation of any jet observable, only the finite terms from each partonic channel will appear. 

As a check, we compute the perturbative corrections to the hadronic i?-ratio, which is obtained by 
simply setting all measurement functions Fj™ 1 ^ to unity, 
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This result is in agreement with the literature [44] , thus providing a further strong check on the correctness 
of all terms computed here. 



6 Summary and Conclusions 

In this paper, we analytically examined the infrared singularity structure of two-jet production in e+e 
annihilation. 
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For this purpose, we developed a subtraction formalism including double real radiation at tree level 
and single real radiation at one loop. In the case of the two-jet production process considered here, the 
subtraction terms coincide with the full tree level four-parton and one-loop three-parton matrix elements. 
The phase space structure of the subtraction extends the NLO dipole subtraction formalism [18] to NNLO. 
In particular, the subtraction terms are integrated (both numerically and analytically) over the inclusive 
phase space. For single-particle subtraction at tree level and one loop, this procedure is already sufficient to 
transfer all infrared singularities from the numerical integration (involving the jet definition) to an analytic 
expression which is canceled against infrared poles from virtual corrections. This procedure has however to 
be extended in the case of double-particle subtraction. Indeed, the inclusive integral of the subtraction term 
will yield infrared singularities from double real emission (which subtract the corresponding terms in the 
full tree level double real emission cross section) as well as infrared singularities from single real emission. 
The latter have to be subtracted from the subtraction term in order to avoid spurious infrared poles in the 
cross section. The infrared singularities from this second level subtraction can in the present calculation be 
identified to correspond to three-jet final states, which are part of the inclusive four-parton cross section. 
They are canceled by explicit infrared poles present in the virtual one-loop single emission subtraction 
term, which is integrated over the inclusive three-parton phase space. In fact, this one-loop single emission 
subtraction term compensates the infrared singularities of the one-loop three-parton contribution to two-jet 
final states, but also yields purely virtual infrared poles in the three-jet region. 

The infrared pole structure of each partonic final state contribution to two-jet final states is computed 
explicitly by analytical integration of the appropriate subtraction terms. To analyze the resulting infrared 
poles and their cancellation among different contributions, we start from Catani's infrared factorization 
formula for two-loop amplitudes [36] and identify individual terms in this formula with corresponding 
terms in the integrals of the tree-level double real emission and one-loop single real emission subtraction 
terms. It turns out that an identification up to the double pole terms is possible, provided a particular 
finite constant (which is irrelevant for the infrared structure of the two-loop virtual corrections) is chosen 
for the infrared singularity operator. Besides these terms canceling with infrared singularities in the two- 
loop virtual corrections, we identify one particular contribution which cancels between double real emission 
and one-loop single real emission. This contribution is found to be proportional to the one-loop correction 
to the soft gluon current [20]. These identifications may help to find a more detailed insight into the 
structure of NNLO infrared singularities of processes with more than two jets in the final state. Once 
this is accomplished, it appears feasible to construct NNLO infrared subtraction terms similar to the ones 
proposed here for processes with higher partonic multiplicity and/or combined emission in initial and final 
state. 
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A Master Integrals 

All loop and phase space integrals in Section ^ were carried out by first reducing them to a small set of 
master integrals. This reduction is based on an iterative solution of integration-by-parts identities [45] and 
is explained in detail in [33,40,41]. For reference, we collected the analytic expressions (which can be found 
in various places in the literature [33,46,47]) for all master integrals in this appendix. 
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A.l Virtual two-loop corrections 



The virtual two-loop vertex master integrals were first derived in [46] in the context of the calculation of 
the two-loop quark form factor [39]. Factoring out a common 



they read 
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A.2 Three-particle phase space integrals of one-loop corrections 

The inclusive three- and four-particle phase space integrals appear in this calculation in the form of dipolc 
and tripolc phase space integrals. We therefore factor out not only the usual normalization factor of 
dimensional regularization, but also the volume of the two-particle phase space into 
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The master integrals appearing as three-particle phase space integrals of one-loop matrix elements are 
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The master integral Vg was first derived in [33] . 



A. 3 Four-particle phase space integrals 

Inclusive four-particle phase space integrals of tree level matrix elements can be expressed as linear combi- 
nation of four master integrals, which were derived in [33]: 
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